FINAL TECHNICAL REPORT 


i 1 --' / ' , s < - • (. -T V 




/// 


^ / 


on 


NASA Grant, NAG 1-1074 


APPLICATION OF OPTICAL DISTRIBUTED SENSING 

AND COMPUTATION 

TO CONTROL OF LARGE SPACE STRUCTURES 


Submitted: May 1992 


A. V. Balakrishnan, Principal Investigator 


(NASA-C*-190203> APPLICATION OF OPTICAL 
OISTs<IdUT?0 SFNSING AND COMPUTATIO 
CONTROL OF LARGE SPACE STRUCTURES f.oal 
Technical Report (Cal i f orni a Un. v . ) 35 p 


CSCL 20F G3/74 


N92-23567 


Unci as 
0086919 



Optical Distributed Sensing for 
Controlling Large Space Structures 

Weijian Zhang 



1 Introduction 

Rapid developments in optical materials and devices in the last two decades 
have made optical sensing a very attractive sensing technique in control en- 
gineering. Compared with the traditional electronic and electro-mechanical 
sensors, the optical sensing and related information processing has the fol- 
lowing advantages: 

1. Interconnections: Optical connections are- much easier to make than 
electronic ones, because optical signals do not affect each other if their 
paths cross. 

2. Immunity to electromagnetic interference (EMI): Electric conductors, 
unlike light beams, when placed too close to each other have EMI 
upon each other, because the signals they carry generate magnetic fields 
which can induce currents in nearby conductors. For example, a prime 
attraction of optical fibre sensors in aerospace is their immunity to 
EMI, including some electronic countermeasures. 

3. Optical measurements are made without making contact, at high speed 
and with high precision. 

4. Reduced size and weight of devices, higher processing capacity for low 
electrical power consumption, are attractive factors to space industry. 
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Although optical fibre sensors have been widely used in industry, their 
applications are limited to the measurements of the following field variations: 

1. strain distribution in large structures such as bridges, dams, and pres- 
sure vessels; 

2. temperature distributions in equipment such as power transformers, 
electrical generators, boilers and high-voltage cables; 

3. electric and magnetic field anomalies in electrical networks. 

In this report a real time holographic sensing technique is introduced and 
its advantages sure investigated from the filtering and control point of view. 
The feature of holographic sensing is its capability to make distributed mea- 
surements of position and velocity of moving objects, such as a vibrating 
flexible space structure. This study is based upon the distributed parameter 
models of linear time-invariant systems, particularly including the linear os- 
cillator equations describing the vibration of large flexible space structures. 
The general conclusion is that application of optical distributed sensors brings 
gain in the situation where Kalman filtering is necessary for state estimation. 
In this case, both transient and steady state filtering error covariance become 
smaller. This in turn results in smaller cost in the LQG problem. 

Instead of measuring field variations like optical fibre sensors, holographic 
sensors are used to measure the position and velocity of objects in motion. 
Rate measurement is made possible by latching two holographic memory 
devices, where the images are retained, at different times. 

The organization of this report is as follows: 

In Section 2, the princple of a real time holographic sensing technique — 
Degenerate Four Wave Mixing (DFWM) — is reviewed briefly. 

In Section 3, we first study noise disturbances in holographic sensing 
and propose white noise as a reasonable formulation of observation noise in 
holographic sensing. 
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Then the results on the gains of optical distributed sensing are stated 
and proved. Finally, by studying the active damping problem of flexible 
structures, we show that it is important to use distributed sensing in order 
to reduce the unwanted energy input caused by filtering errors. 

Conclusions are in Section 4. 


2 The Principles of Real Time Holography 


The principle of holography 

(1) Recording of hologram. For simplicity, let us consider only the one dimen- 
sional case. The arrangement of hologram and coordinate is as illustrated in 
Figure 1. 

Obviously, when |x| <C |z,|, i — 1,2, we have 

r(x, *.) = (*? + x y» = M + \ + o ((| )’) . 

The reference wave W rt j(x, z), is a plane wave and is coherent with the 
object wave W 0 f,j(x, z). In phasor form, we have 

W re j(x,z ) = I< exp [ifc (xcos (f — + ycos f + zcosfl)] 

= K exp[i/c(:r sin 6 + z cos «)] 

where k = 27t/A is wave number and A is wave length. 

For simplicity, we consider an object consisting of only two points with 
spacing \z x — z 2 [. Then the object beam is the inteference of two spherical 
waves produced by the two points. Therefore the object wave in phasor form 
is given by 


[ 2—0 


e ikr i z 1 1 *) 
a r(z lt x) 


+ a 


e t~fcr<» 2 ,x) 

r(z 2 t x) 


« a c^i e .fcrV2|z 1 | +a £i^i e .>xV2b 2 | 


= aAuP*™*' I) + aA 2 e‘**V(2M) 
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"''Reference wave (a coherent plane wave): j{ e ,k ( zco * 8 + xt,n6 ) 


Figure 1: The Schematic for Recording a Hologram 
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for |x| < |z 2 |. 

If we assume the exposure E is in the linear range of the transmittance-exposure 
relation [3], then what is recorded on the hologram (x-axis) is the intensity of the 
interference pattern between the object wave W 0 bj(x,z) and the reference wave 
W re j{x,z) at z = 0. 

I(x) = \W ref (x,0) + W oij (x,0)\ 2 

= | Ke ikx ' in$ + aA 1 e* fc * a/(2| * ,|) + aA 2 e''* x:,/(2|x:,,) | 2 
= I< 7 + a/fe , ' fcx, ^*(^e“ , ' fcx3/(2| ' ,l) + A 2 e" il:x2/(2|xil) ) 

+a/fe- ,fcx " ntf (v4 1 e ,fcxJ ^ (2| * l|) + A 2 e’* xJ/ < 2|,j|) ) 

+0(a 2 ) (1) 


for K a. 

(2 )Reconstruction of the recorded holographic image. Those two cross product 
terms in (1) play crucial role in holography. After developing the film, put it in the 
reconstruction system as illustrated in Figure 2. 

Let us consider the contribution of the first cross product term. The spatial 
domain solution of the diffraction problem is expressed in the Fresnel-Kirchhoff 
integral as follows: If a plane wave of amplitude a 2 , traveling in the direction of the 
positive z axis, is incident on an object with amplitude transmittance £(xx,yi) in 
the plane normal to the z axis at z = 0, the light complex amplitude a 2 (z 2 ,y 2 ,d) 
in the plane z — d is 


< 12 ( 3 : 2 , 2/2 ,d) = 


exp{— t(27r/A)[d 2 -f (x 2 - xx) 2 + (y 2 - r/x) 2 ] 1/2 } 
X [d 7 + (x 2 - xx) 2 + (y 2 - yi) 2 ] 1/2 


dxidyi 


The geometry is illustrated in Figure 3. The angle 9 is formed by the positive z axis 
and the straight line connecting the points (xx,yx,0) and (x 2 ,y 2 ,d). Cos0 is called 
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Figure 2: The Schematic for Reconstructing the Image 
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the obliquity factor , which, in what follows, is considered to be close to 1 since the 
angle 0 is generally close to 0 when d is large. 

By Fresnel- Kirchhoff integral [3, pl06], and neglecting the multiplicative con- 
stant, we obtain the diffraction pattern of the read beam upon the recorded holo- 
gram as 


9(0 


e ifc(*+< 3 /(2*)) J ( fil e-** 3 /(2|*l|) + J 4 2e -ife* a /(2|* 3 |)) 

x exp[iA;x 2 /(2z) — ikx/z(£ — zsin0)]dx 


Therefore, at position z = |z 2 |, the first term in the above integral evaluates 
to 6(£ — |zj|sin0), while the second term becomes — |z 2 |sin0) at z = |z 2 |. In 
other words, the two point object is reconstructed at positions (|zi|, |zx| sin0) and 
(|z 2 |,|z 2 |sin 0). 

A Real-time Holography Technology - DFWM 

In present days, real-time holography can be accomplished in photorefractive 
materials, a technique called degenerate four wave mixing (DFWM) [4]. Photore- 
fractive materials include the following type of crystals: lithium niobate (LiNbOa), 
potassium niobate (KNbOa), barium titanate (BaTiOa), strontium barium niobate 
(SBN) and bismuth silicon oxide (Bii 2 Si0 2 o). Photorefractive crystals provide a 
unique way of recording light intensity which allows multiple beams to be mixed. 

To understand the operation of DFWM, let us first recall that the refractive 
index n, or index of refraction, of a medium is defined to be c/v , where c is the speed 
of light in free space and v the speed of light in the medium [7]. For photorefractive 
crystals, their refractive index is nonlinearly dependent upon the local intensity I 
of light incident upon it. The total refractive index is given by [9] 

n = n 0 + n 2 / p (2) 


8 




Figure 3: Geometry for Fresnel-Kirchhoff integral 
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where no,n 2 ,p axe constants and it is common that p ^ 1. 

When light is incident upon a photorefractive crystal, the refractive index 
within the crystal changes with the local light intensity, by (2). In this way, 
the distribution of light intensity over the volume of the crystal is recorded 
as a change in refractive index. 

When performing DFWM using photorefractive crystal, as in conven- 
tional holography, both the object beam and the reference beam, which axe 
coherent, are incident upon a photorefractive crystal. A1 the same time, a 
third beam, called read beam, is also incident upon the crystal, see Fig- 
ure 4. The interference pattern generated by the interference of these three 
beams produces a spatially varying refractive index within the crystal. The 
variation in refractive index within the crystal is similar to a recorded holo- 
gram, and causes diffraction of the beams. If the read beam is identical but 
phase conjugate to the reference beam, then the diffracted beam is the phase 
conjuate of the object beam. 

DFWM in a photorefractive crystal is similar to conventional holography 
where the hologram is illuminated with the phase conjugate of the reference 
beam, producing a phase conjugate of the object. The difference between 
DFWM and conventional holography is that in the DFWM process, the 
hologram is made and read simultaneously, avoiding the chemical process of 
developing the hologram. Because of the fast response times of many pho- 
torefractive crystals [5], the holography can be made real-time using DFWM 
in a photorefractive crystal. 

Currently, there are two major factors that limit widespread application 
of photorefractive crystals: 

1. Some of the most promising crystals (e.g., BaTi0 3 , SBN and KNb0 3 ) 
are not widely available in large samples with high optical quality. 
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2. The crystals that are available do not perform well enough in all re- 
spects. For instance, in order to demonstrate high speed, low write 
energy, long memory or large gain, it is necessary to use several differ- 
ent types of materials. Reference [11] gives the performance evaluations 
of some photorefractive crystals. 

3 Results 

Noise disturbances in holographic sensing 

Noise disturbances are always present in the processes of making a holo- 
gram and then reading the hologram. Therefore we have to consider sensor 
noise in optical holographic sensing. Generally, the noise is attributed the 
following sources: 

1. Random scattering of both signal and reference beams during exposure 
due to the granularity of the recording material. 

2. Random scattering of the reading beam and the reconstructed wave 
due to the granularity of the material of the hologram. 

3. Spatial modulations of the reference and reading beams. 

4. Inhomogeneities and surface deformations of the recording material. 

5. Nonlinear recording of the signal wave. 

Sources 1, 2 and 3 are usually not significant except in certain special 
cases such as making multiple recording of many holograms in the same 
material, or using spatially modulated reference beams. Source 4 is related 
to the preparation of the recording material. Next we study the adverse 
effects of source 5, the nonlinearity in holographic recording. 
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Figure 4: The scheme of Degenerate Four Wave Mixing (DFWM) 
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The transmittance t of a hologram made is a function of the exposure E, 
which is defined by 

E = kiT e I 

where ki > 0 is a constant, r e is the exposure time. In the previous analysis 
we have assumed that the relation between t and E is linear. However, 
t — E curve is always nonlinear to some degree. The following form of cubic 
approximation 

t = Co + ciE + c 2 E 2 + C3E 3 (3) 

fits experimental curves quite well over a considerable range for plane ab- 
sorption hologram. 

For the usual range of exposure, nonlinear effects produced by the quadratic 
term are the most important. Its nonlinear effects are particularly well illus- 
trated for a two point object. If we adopt the reconstruction arrangement as 
in Figure 5, then the wave at z = d (on the hologram) is given by 

« = W obj (x2,y 2 ,d) 

= exp (-ikyjxl + (y 2 - W) 2 + d 2 ) + exp [~iky/x\ + (y 2 - & 2 ) 2 + ~d 2 ) 

~ exp [-ik / (2 d)(xl + yf - 2y 2 6!)] + exp [~ik/(2d){xl + y\ - 2y 2 6 2 )] 

= e '^ 1 + e '* 2 

For simplicity of analysis, let us assume that the reference wave is unmod- 
ulated, i.e., rr* = constant. Then the quadratic term in the transmittance 
of the recorded hologram is given by 

c 2 E 2 = c 2 k?T e 2 ( |a + r| 2 ) 2 
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Figure 5: The generation of false images of a two point object by a nonlinearly 
recorded hologram (h = £>j — 6 2 ). 
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= c 2 fcjr 2 (aa* + rr* 4- ar m 4- a*r) 2 

= c 2 fc 2 r 2 [(aa*) 2 + (ar*) 2 + (a*r) 2 + 2 aa*ar* + 2aa*a*r] 
+L(aa*,ar*, a*r) 

where Z/(>, •, •) is a linear function of its three variables. 

When the nonlinearly recorded hologram is illuminated by the original reference 
wave r, the fourth term in the above equation produces the following diffracted 
wavefront 

r(2aa“ar“ ) = 2|r| 2 (a 2 a*) 

= 2 | 4 - + e~* 3 ] 

= 2|r| 2 [3e‘ >1 + 3e*^ + c *Wi-*a) + e «'(2*»-*i)] 

= 6|r| 2 (e*^ 1 + e’* J ) 4- 2|r| 2 [exp(-tfc/(2d)(x 2 4- y\ - 2y 2 (bi + h ))) 

+ exp(-ik/(2d)(xl + y\ - 2y 2 (6 2 - h))) ] 

where h = b\ — b 2 > 0. 

It can be seen that the third term in the above represents a spherical wave 
diverging from a false virtual image of a point source Pa(0, b\ + h). And the fourth 
term in the above represents a spherical wave diverging from another false vir- 
tual image of a point source /^(O, b? — h). Therefore, two false virtual images are 
produced by the effect of nonlinearity in t — E relation. 

The cubic nonlinearity term in (3) can be similarly analyzed and additional false 
images are introduced in and around the true image. 

Therefore, there are numerous sources of noise in holographic sensing and they 
are in and around the true image. We propose the white noise formulation to 
represent the observation noise. 


15 



The gain of optical distributed sensing 

Let us consider the following stochastic system 

i(t) = Ax(<) + Bu(t) + Fu>(t) 
y(t) = Cx(t) + Gu>(t) 

where GG‘ = 7, FG m = 0 (noise independent), and B is Hilbert-Schmidt, and {*?(•) 
is a Gaussian white noise process. 

In the associated LQG problem, the cost functional is given by 

J(u) = lim f Tr.QE[x(t,u))x‘(t,u)]Q‘dt 

T—oo i Jq 

+ / Tr.£;[u(t,u;)u*(t,u;)]dt} (4) 

1 Jo 

where Q is a Hilbert-Schmidt operator. 

By [1], the optimal feedback control for the LQG problem is given by 

u*(t) = —B m P c x a (t) 

where P c satisfies the following SSRE 

A‘P C + P e A + Q‘Q - P C BB'P C = 0 (5) 

and x a (t) is defined by the asymptotic Kalman filter 

i«(<) = (A-BB'P c )x a (t) + PjC-(y-x a (t)) 

= (A - BB‘P C - PjC~C)x a (t) + PjC‘y(t) 

£ o (0) = 0 

where Pj solves the following SSRE 

APj + PjA m + FF‘ - PjCCPj = 0 (6) 
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And the minimal final cost is given by 


J(u) = Tr.QPjQ * + Tr.CPjP c PfC m (7) 

Sufficient conditions on the existence of P c , P/ such that ( A — BB m P e ) and 
{A — P S C'C) are strongly stable 

1. For our particular 

a= {-a o) ; B= (b) 

by [2], if (A, B) is controllable and 

||Qu>|| 2 < M||0*u>|| 2 Vt oeV(VA)xH 

and also (.4.*, Q m ) is controllable, then (5) has unique solution such that (A — 
BB*V C ) is strongly stable. 

2. By [2], if ( A*,C *) and (.4, !F) are controllable and 

ll^*HI 2 < M||Ctu || 2 


then 

Hm Vj(t) = Vj (in strong sense); 

and Vj satisfies (6) and is the unique solution such that both 
(A — V/C*C) and (.4 — VfC m C)* are strongly stable. 

Zabczyk [12] obtained more demanding sufficient conditions which require the 
exponential stabilizability of (A, B) } (,4,.F), (,4*,C*), (^4*,Q*). Unfortunately all 
these conditions are not satisfied in view of Gibson’s result [6] and the fact that B 
and Q are Hilbert-Schmidt. 



From the above well established results, we know that no matter what 
kind of choice of the observation operator C we make, the optimal control 
law is given by 

«*(<) = -B*V c w a (t) 

where "P c is independent of the observation operator. Therefore the gain 
of distributed sensing will only happen in filtering, in the sense that V/ = 
lim t - t0 o'Pf(t) will be smaller. Furthermore, as a consequence of “better” 
recovery of state information by Kalman filtering through distributed sensing, 
the optimal feedback control should give a smaller final cost. 

Next we first consider the following two SSRE’s in Kalman filtering. 

AP + PA* + FF* - PC*CP = 0 (8) 

AP d + P d A* + FF* - P d C* d C d P d = 0 (9) 

where C : H — *• IR n , C d : H — ► H 0 are linear and bounded. 

We are going to show 


Proposition 1 Under the assumptions 


CC d C d = CC*C 
CC* > 0 


there holds 


P d <P 


provided that both P and P d exist. 


( 10 ) 

( 11 ) 


Before the proof of the proposition, we notice that the assumption (10) 
has the clear physical interpretation that the information yielded by the 
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observation operator C is just a subset of that of G d - 

Proof: Consider the following two LQR problems on the infinite horizon, 


f x(t) = A"x(t) + C*u(t) 
( x(0) = x 


with cost functional 


roo 

J(u;x)= / ||F*x(0ll 2 + IK0l| 2 Kndf 

JO 

f x(t) = A*x(f) + G d u(t) 
x(0) = x 


with cost functional 


M“;*) = II^WII 2 + IK<)llk <«• 

For P and Pd satisfying (8) and (9) respectively, we know that 
\Px,x] = min J(u:x ) 

uE£/ 2 [(0 t oo),H n ] 


[P d x y x)= min J^ujx). (15) 

u6Zy 2 l(0 t OO),tt O ] 

Next, in the second LQR problem (13) we perform optimization with 
respect to only a special class of control. 


Since 


u = CdC'iCC'y'u with u £ L 2 [(0, oo), JR n . 


C' d u = C2C d C'{CC*)-'u 


= ( CC' d CdY{CC')-'u 

= (cc*cy(cc*)- l u 


= C’u 
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and 


IMIk = [cc;c J c-(cc-)-'«,(cc-)- , u] 
= INI’*. 


then the second LQR problem becomes 

f x{t) = A-x{t) + Qu(<) = A*x(f) + C'u{t) 

y x(o) = x 

with cost functional 

fOO „ _ 

UCtC-icc-y'u-,!) = / ||f x(()|| 2 + iiwtccr'ailk* 

Jo 

= / o °° iif *(on* + iwoift** 

= J(u-,x) (16) 


Therefore, in the light of (1 4)-( 1 6) , we obtain 


[PdX,x] = 


< 


min JJu;x) 

ueL s ((0,oo).«ol 

min J d (CdC*u-,x) 

ueJ^Ko.ooj.w] 

min J(u:x) 

u€L j [( 0 ,oo),R"] 


[Pi, x] 


□ 

Remark: Notice that in the proof of Proposition i, we only need the assumptions 
(10) and (11), and C d is not necessarily a distributed sensing operator. In addition, 
the requirement that Range(C) = 2R n - a finite dimensional space, is superficial. 
Therefore, we have proved the conclusion: 


CiC;c 7 = & CiC; > o => p }a < p } , i 


i.e., more observation implies smaller steady state filtering error. 
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It is also interesting to apply the above conclusion to finite dimensional 
Kalman filtering. 

Next, we study the influence of distributed sensing to the value of final 
cost. 


Proposition 2 In addition to the assumptions (10) and (11), suppose (A, B) 
and ( A *, Q *) are exponentially stabilizable. Then, 


min 

u(- I a>)eL J ((0,oo),H„] 


Jd{u) 


min 

« ( * »<^) € L 2 [(0 ,oo) , H u ] 


J(u) = Tr. B'P c (P d - Pj)P c B < 0 


(17) 

where «7(u) and J d (u), as given in (4), are the quadratic cost functionals 
corresponding to the systems with observation operators C , C d , respectively. 


Here we impose the much stronger conditions of exponential stabilizability 
of (A,B) and (i4*,<5*) as a sufficiency condition such that ( A — BB* P c ) is 
exponentially stable. And hence 

J r oo 

1 S;(t)Q*QS(t)x dt xeH 
0 

which implies P c is nuclear, where S(t) and S c (t) are the Co-semigroups 
generated by A and (A — BB*P C ), respectively. 

In general, P c need not be nuclear without the above assumptions. For 
example, in the usual notations, let 



where B : H u — ► H is Hilbert-Schmidt. Then, 


A + A" = 0 => V c = I 

which is not even compact. Obviously, in this particular case, ( A,B ) and 
are not exponentially stabilizable in view of Gibson’s results [6], 
because B and Q* are compact. 


21 



Proof: Recall that the minimum cost of (4) is given by (7) and Pf and P c satisfy 
(6) and (5), respectively. 

Let {e n } be an orthonormal basis in H . Then the minimal cost can be rewritten 
as 


J(tx‘) = Tr.QPjQ’ + Tr.CPjP c PfC m 
= Tr.PjQ'Q + Tr.PfC'CPfPc 

= ^{[Q'Qt^PjcA + lPiC'CPjPctn,^]) 


n= 1 
oo 


= £{[P/(P c BB*P c - A* P c - P c A)e n ,e n ] 

n= 1 

+[(AP f + PjA * + F F*)P c e n , e n ]} 

OO 

= E {[P/PcBB'P c e n ,e n ] + [FF‘P c e n ,e n ]} 

n=l 

= Tr.PjP c BB" P c + Tr.FF'P c 
= Tr.B m P c PjP c B + Tr.F*P c F 


Similarly, we have 


Jd{vT) = Tr.B*P c P d P c B + Tr.F*P c F 

Therefore, subtraction of the two gives (17). □ 

The conclusions of the above propositions can be extended to the finite time 
horizon LQG problem: 

£(f) = Ax(t) + Bu(t) + Fu>(t) 0 < t < T 
y(t ) = Cx(t) + Gu(t) 

Ex(0)x(0y = A (18) 
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where GG * = /, FG * = 0, and B : J/ u — > H is Hilbert-Schmidt. The cost 
functional is given by 

rT [T 

J(u) = I Tr.E[(Qx(t,w)(Qx(t, + / Tr.E[u(t,u})u(t, u>)*]df 
Jo Jo ' 

where Q : H — ► H is again Hilbert-Schmidt. 

Let us consider the differential Riccati equation corresponding to the Kalman 

filtering 

P f (t) = APj(t) + P f (t)A- + FF- -Pf(t)C m CP f (t) 

Pfi 0) = A (19) 

and the differential Riccati equation corresponding to the optimal feedback gain 

Pc(t) = -A‘P c {t)-P c (t)A-Q’Q + P c {t)BB 9 P c (t) 

P C {T) = 0 (20) 

In what follows, we use Pj(-) to denote the solution of (19) with C being replaced 
by Cd , and Jd{u) is similarly defined. Then we have 

Proposition 3 Under the assumptions 

CCjCd = cc m c 
CC m > 0 

there hold 

Pd{t) < P/( t) for 0 < t < T 


min 

U6L»[(0,T);// U ) 


Jd(u) 


min J(u) 

ueL*[(0,T)-,H u ) 


= [ T Tr.B-P c (t)[P d (t) - Pj{t))P c {t)Bdt < 0 

Jo 
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Proof: First of all, By [1] we know Pd(’) and F c (-) all exist and are unique. 

Let K(t) = P/(T — t ) and K d (t) = P d (T — t), then K(t) satisfies 

K{t) = -AK(t) — K(t)A* — FF* + K(t)C*CI<(t) 

K(T) = A 

while K d (t) satisfies the same equation with C being replaced by C d . 

Consider the following finite horizon LQR problems 

x(£) = A"x(t) + C*u(t) 

x(t 0 ) = x 0 < to < T 

J(W,u,x) = (A*(T),*(r)] + l T [FF-x(t),x(t)] + ||i(i)H 2 <<( (21) 

J to 

and 


x(t) = A m x(t) + CJu(t) 
x(t 0 ) = x 0 < to < T 

Jd{to\ u,x) = [Ax(r),x(r)]+ f T [FF'x(t),x(t)] + \\u{t)\\ 2 dt (22) 

Jto 

It is well known that [1] 

min J(t 0 :u,x ) = [KUq)x,x\ = [Pt(T — to)x,x] 

u€lPI(0,T) t H"] v ' *" 1 J 

mjn „ , ^(* 0 ; x) = [j P d (T - t Q )x , x] 

u€L 2 [(0 t T),H o ] 

Then, similarly as in the proof of Proposition 1, by taking 

u = C d C m {CC')~ x u 


we can establish 


[Pd(T - £ 0 )a:,a:] 


< 


min 

u€I,»[(O.T),H.J 


Jdit 0 \ u,x) 


min J d {t 0 ,C d C'{CC‘)~ l u,x) 
[Pj{T - t 0 )x,x] 
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For the second half of the conclusion, we can show in a similar manner as in the 
proof of Proposition 2 , that 


min J(u) 
ueL*[(0,T),Hu) 

= [ T Tr.B'P c (t)Pf(t)P c (t)B + Tr.F P c (t)Fdt + Tr.\P c (0) 
Jo 


from which the result follows by subtraction, noticing that P c (-) is independent of 
the observation operator C or Cj. □ 

Next, let us consider the effect of distributed sensing on active damping of 
large flexible space structures. Since we are concerned with the design of stablizing 
control law instead of numerical solution of the structure response, we use the 
following continuuam model in order to avoid losing the insight of the infinite 
dimensional nature of flexible structures, 


x(f) + 7 I>(x(<), *(*)) + Ax(t) = Bu + FN(t) (23) 

where A : T>(A) — ► H is self-adjoint, nonnegative, has compact resolvent and 

dense domain. We further assume that 0 belongs to p{A), the resolvent 
set of A , i.e., there is no rigid body mode. 

B : M n ► H is linear and bounded. 

F : H n — ► H is a Hilbert-Schmidet operator. 

£>(x, x) is the (generally nonlinear) damping term, such that [£>(x,x),x] > 0. 
N(t) represents a white noise in the Hilbert space H n . 

Since the damping coefficient 7 > 0 is typically very small in flexible space 
structures, we set 7 = 0 in the following studies for the sake of simplicity. 

As usual, we may define 
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so that (23) can be rewritten as 

- = AY(t) + Bu + TN{t) where Y = ^ (24) 

with T>(A) = T>(A) ® H and equipped with the energy norm 

||K|||=[/l ,/ V>l ,/J x| + [x,i] 

on the Hilbert space He = T>(A X ^ 2 ) ® H. 

From practical point view, only finite number of controllers can be present 
and they are of saturation type due to limited control forces and moments [10]. 
Therefore, 

ti 2 

where £/,* are certain positive constants derived from maximum force/moments avail- 
able. 

The total energy of the flexible structure governed by the continuum model 
(23) is defined by 

E(t) = 1/2||F(0IIe = l/2([Aar(0.x(01 + [i(<).*(0D 

Since the energy decay rate is given by 

^ = [u,B-x(!)]+[™(f),i(()] 

then obviously the optimal feedback control for maximum energy decay rate is given 
by 

Uj{t) = -UjS\gn[^p-] j = 1, 2, • • • , n 
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where 


b(t) = 


l m \ 
MO 


= B m x{t) 


\ MO ) 

Therefore the control law only requires the rate information at the locations of the 
controllers b(t) — B*x{t). In other words, the rate measurements made by colocated 
sensors 


y(t) = B-x(t) 


are necessary and also sufficient for active damping purpose. More observation than 
that made by colocated sensors does not make any difference in terms of energy 
decay. Therefore, in this particular case, optical distributed sensing does not bring 
any gain. 

However, if we take sensor noise into consideration, then the situation will be 
quite different. 

Next, let us consider the case in which the observation of structure deflection 
rate is corrupted by additive sensor noise, i.e. 


v(t) = Cx{t) + GN(t) 

where FG* = 0 (observation noise and state noise are independent), and C could 
be colocated sensor B* or optical distributed sensor Cd for later comparison. 

For simplicity of analysis, let us assume that linear feedback is allowed, which is 
reasonable for small deflection magnitude. Then a reasonable choice for stablizing 
feedback control would be 


u(t) = 


Then the closed loop system dynamics becomes 

IXiH = AY{t) — BB‘Y(t) + FN(t) 
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= (A - BB')Y(t) + BB m (Y(t) - Y(t)) + FN{t) 

We therefore see that due to filtering error, we thereby introduce a random input 
to the system which may excite higher order modes. And in this case, it is easy to 
verify that the energy decay rate becomes 


^ = — IIB-itOII’ + l|fl'(i (0 - i( 0 )ll* 


(25) 


The second term in (25) represents a positive energy input, which discounts the 
active damping provided by the controllers. The third and the fourth terms can 
be considered as zero mean disturbance. Therefore, we are facing two options to 
further enhance the structure stability: 

1. Increase the active damping by putting more controllers and increasing the 
force/moment upper bounds, i.e., to increase the stability margin of ( A — BB ) 
by changing B . But this option generally results in heavier weight and bigger 
cost. 

2. Decrease the positive energy input 

||B-(±(t) - x(l))H 2 

which is caused by filtering error, i.e. to decrease P 4 , the error covariance 
operator of x(t) — x(<), by employing optical distributed sensor. 

Let 

v= {p 3 p.) 

be the steady state Kalman filtering error covariance operator. It is generally a 
formidable problem to solve the corresponding SSRE 


A V + VA * + FF - VC'CV = 0 


(26) 
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to find V. 

Next, a very wild attempt is made to find an “approximate” solution of (26). 
First of all, since we are interested in a self-adjoint solution, it is necessary to 
find the expression of V m with respect to the PLe inner product. It is easy to verify 
that 

_ / A~ x PfA a-'p ; \ 

V p l A p < ) 

Then "P is self-adjoint if and only if Pj, j = 1,2, 3, 4 satisfy the following equations 


APi = (AP l y (27) 

Pz = (APiT (28) 

Pa = P; (29) 

As above, suppose only rate measurement is made, i.e. C = (0,C), then direct 

computation using (26) gives 

p 2 + p ; - PiC'cp; = o (30) 

P*- Pi- P 2 C‘CP 4 = 0 (31) 

APi + P;A- FF* + P 4 C m CP 4 = 0 (32) 


Therefore, to solve (26) for a self-adjoint solution is equivalent to solving equa- 
tions (27) - (32). For this purpose, we first eliminate Pi and P 2 from (27) - (32) by 
noticing 

Pi = A~'P; (33) 

P\ = P 4 - A~ l P;C‘CP 4 (34) 

from (28) and (31). 

Left- and right-multiplication of A to (30) gives 

ap 3 + p;a - p;c*cp 3 = o ( 35 ) 
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And (32) can be rewritten as 


P 3 + P; - FF m + P 4 C~CP 4 = 0 (36) 

Substituting (34) into (27) gives 

ap 4 + P 4 C-CP 3 = p 4 a + p;c*cp 4 (37) 

Therefore, the problem becomes solving Pz and P 4 from (35) - (37) and (29). 
Once P 3 and P 4 are found, P x and P 2 are immediately given by (33) and (34). 

Now we attempt to find an “approximate” solution of the form P 3 — 0, in which 
case, P 2 = 0 and Pi = P 4 . 

By the assumption on A , we know there is a sequence of normalized modes <f> n s 
and corresponding frequencies u> n ’s such that 

A<f> n = u 2 n <f> n n = 1,2, •• • 

with lim n _oo — 00 and { <f> n , n = 1, 2 , • • • } being an orthonormal basis in H. 

Since P 3 = 0, (35) is trivially satisfied. We only need to determine P 4 from (36) 
and (37), which become 

FF m = P 4 C'CP 4 (38) 

AP 4 = P 4 A (39) 


and (29). 

If we assume iv n ’s are distinct, then it is easy to show that 

AP 4 = P 4 A and P 4 = P 4 

if and only if there is a sequence of nonnegative numbers { a n , n = 1, 2, • • • } such 
that 

P 4 <j> n = n=l,2, ••• 


) 
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) 


Now a n ’s can be determined from (38), i.e. 

{FF'4> n ,<l> n \ = [P 4 C*CP 4 <t> n ,<f> n ) = a^||C7^ n || 3 

which immediately gives 

llfV.ll „ _ 1 o 


a„ — 


\\C*n\\ 

Therefore, we have found an “approximate” solution V of (26) given by 


■(v ;.) 


where P 4 is a self-adjoint, nonegative operator defined by 

p , IIF-A.il, . , 

i^7T T T » n — 1|A" - 


with domain 


iicA .11 

V(P t ) = {* € H | £(1£M) J [*, Kf < °° } 


(40) 


S' IICA.II 

This is an “approximate” solution in the sense that it satisfies only equations 
(35), (37) and (29), but not (36). In fact, if we consider the left hand side of (36) as 
an infinite dimensional matrix, then the above solution only makes all the diagonal 
elements zero, while the off diagonal elements are not necessarily zero. 

It is easy to see that a sufficient condition for (40) to be an exact solution of 
(26) is 


[FF’<t> n , <j> m ] = [P 4 C*CP 4 <t>n, K] = a n a m [C'C<t > n , <£ m ] Vn # m 


or, equivalently, 

IfVnyf^v) jCAnCAj y , 
IIF-A.il IIF-A.il IICA.II IICA.II r 

which is a very strong requirement. 


) 
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4 Conclusions 


In summary, optical distributed sensing brings gain in the situation where Kalman 
filtering is necessary for state estimation. In this case, both the transient and 
the stationary Kalman filtering error covariance become smaller. Due to smaller 
filtering error, the same control law results in smaller cost in LQG problem. These 
conclusions;, of course, also apply to finite dimensional Kalman filtering and LQG 
problems as special cases. 

By considering the stabilization problem of flexible space structures, we realize 
that it is important to use distributed sensors in order to reduce the positive energy 
input caused by filtering error. 
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